Checkerboard local density of states in striped domains pinned by vortices 
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We discuss recent elastic neutron scattering and scanning tunneling experiments on High-Tc 
cuprates exposed to an applied magnetic field. Antiferromagnetic vortex cores operating as pinning 
centers for surrounding stripes is qualitatively consistent with the neutron data provided the stripes 
have the antiphase modulation. Within a Green's function formalism we study the low energy elec- 
tronic structure around the vortices and find that besides the dispersive quantum interference there 
exists a non-dispersive checkerboard interference pattern consistent with recent scanning tunneling 
measurements. Thus both experiments can be explained from the physics of a single Cu02 plane. 
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The competing orders in the High-Tc cuprates remain 
a strong candidate for explaining some of the unusual 
features of these doped Mott insulators^"^. The compe- 
tition between superconducting order and antiferromag- 
netic order has recently attracted a large amount of both 
experimental and theoretical attention. In particular, ex- 
periments in the mixed state have revealed an interesting 
coexistence of these order parameters. 
Elastic neutron scattering results on La2-xSr2:Cu02 
(x=0.10) have shown that the intensity of the incom- 
mensurate peaks in the superconducting phase is consid- 
erably increased when a large magnetic field is applied 
perpendicular to the Cu02 planes^. This enhanced inten- 
sity corresponds to a spin density periodicity of eight lat- 
tice constants Sag extending far outside the vortex cores. 
Similar results have been obtained for the related ma- 
terial La2Cu04+j,^. Nuclear magnetic resonance (NMR) 
experiments have shown evidence for antiferromagnetism 
in and around the vortex cores of near-optimally doped 
Tl2Ba2Cu06+5*. Furthermore, muon spin rotation mea- 
surements from the mixed state of YBa2Cu306.5o find 
static antiferromagnetism in the cores^. Consistent with 
these findings scanning tunneling microscopy (STM) 
measurements performed on the surface of YBa2Cu307_5 
and Bi2Sr2CaCu208+a:^°'^^liave revealed very low DOS 
inside the vortex cores^^"^*. Thus, there is increasing 
evidence for antiferromagnetic correlations in the vortex 
cores of the under- and optimally-doped regime of the 
hole doped cuprates. More recent STM measurements 
of slightly overdoped Bi2Sr2CaCu2 08+a: have shown a 
checkerboard halo of the local density of states (LDOS) 
around the vortex cores^^. This LDOS modulation ob- 
served at low energy < 12 meV was found to have half 
the period of the spin density wave (SDW) observed by 
neutron scattering (i.e. four lattice constants 4ao), and 
to be oriented along the crystal axes of the Cu-0 plane. 
The neutron experiments have been analysed within phe- 
nomenological models that assume close proximity to a 
quantum phase transition between ordinary supercon- 
ductivity and a phase with antiferromagnetism or a phase 
where superconductivity coexists with SDW and charge 



density wave (CDW) order^'16"^^. In these models the 
suppression of the superconducting order inside the vor- 
tex cores allows the competing magnetic order to arise. 
Demler et ed}'^ found that around the vortices the circu- 
lating supercurrents can similarly weaken the supercon- 
ductivity and induce a SDW. 

The field-induced checkerboard LDOS pattern in the 
mixed state has been recently considered within the 
framework of several models^*"^^. In this paper we add 
to the discussion by calculating the LDOS in regions 
where a d-wave superconductor has been perturbed by 
induced magnetism. First, however, we note that a 
checkerboard spin modulation is inconsistent with the 
elastic neutron scattering experiments by Lake et al.^ on 
La2-2;Sr2;Cu02 (x=0.10). For example, assuming that 
the checkerboard CDW is intrinsic to the Cu-0 planes 
where it gives rise to a static SDW checkerboard pat- 
tern (Fig. la), the expected neutron diffraction pattern 
is shown in Fig. lb. 
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FIG. 1. a) Real space picture of the spin structure in a 
checkerboard spin geometry. Black (white) represent spin up 
(down) while gray reveals the superconducting background. 
In order to simulate the induced incommensurability each is- 
land of antiferromagnetic spins is out of phase with its nearest 
neighbor, b) Fourier spectrum of the spin checkerboard struc- 
ture shown in a). 

As is evident there is a 45 degree rotation of the four 
main incommensurate peaks and a plaid pattern of the 
higher harmonics. The rotated incommensurability (with 
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the correct absence of an increased signal at (tt, tt)) shows 
that this spin structure does not apply to LSCO for dop- 
ing levels close to x = 0.10. It is interesting to note 
that a rotation of the incommensurable peaks at low 
dopings {x < 0.055, close the insulator-superconductor 
phase transition) has been observed in LSCO^^. How- 
ever, there is no simple way to create an antiphase spin 
geometry without frustrating the spins at low dopings 
where droplets of charge in an antiferromagnetic back- 
ground is the expected situation^^. However, this might 
be possible in the highly overdoped regime where the 
droplets have been inverted to separate magnetic islands. 
In that case a 45 degree rotation of the incommensurable 
peaks would be consistent with a checkerboard spin pat- 
tern. In this light it would be very interesting to perform 
an experiment similar to Lake et al.^ on highly overdoped 
LSCO. In the case of a connected antiferromagnetic back- 
ground one would also expect a large weight at (7r,7r). 
The physical picture we have in mind is presented in 
figure 2a. In this real space picture an antiferromagnetic 
core (center) has pinned a number of surrounding stripes. 




FIG. 2. a) The idealized version of a real space spin con- 
figuration consistent with our physical picture, b) Fourier 
spectrum of the spin density order from a). Almost all the 
induced weight is located in the four incommensurable peaks. 



This pinning effect of SDW by magnetic vortex cores is 
a well-known effect from numerical studies^ ^. 
Both experimentally^* and theoretically^'^'^'^^ we expect 
an antiphase modulation of the induced antiferromag- 
netic ring domains. Indeed as seen in Figure 2b the re- 
lated diffraction pattern is qualitatively consistent with 
measurements by Lake et aJ.® of enhanced intensity at 
the incommensurate points. 

Without an applied magnetic field, only disorder can pro- 
duce a similar pinning effect of the fluctuating stripes^''. 
In addition to the creation of more pinning centers when 
applying a magnetic field, the single site impurities are 
expected to pin much weaker than the large "impurities" 
created by the flux lines. This is qualitatively consistent 
with the measurements by Lake et al.^ of the temperature 
dependence of the increased magnetic signal for different 
magnetic field strengths. 

This leads to the question of the electronic structure 

around extended magnetic perturbations in d-wave su- 
perconductors. The many experiments indicating coex- 
istence of d-wave superconductivity and antiferromag- 



netism mentioned above motivate studies of simple mod- 
els that enable one to calculate the LDOS in such regions. 
The model Hamiltonian defined on a 2D lattice is given 
by 

H = — ^ ^ tnm^na^rna- ~ M ^ ^ ^na^na (1) 

{n,7n)(7 na 

^-*=^M„(cl^c„T-ci^c„i) (2) 

n 

where ^ creates an electron with spin a at site n and ji 
is the chemical potential. The staggering is included in 
M„ = (— 1)"M. The strength of the antiferromagnetic 
and superconducting coupling is given by M and A, re- 
spectively. 

The Hamiltonian + is a simple mean-field lat- 
tice model to describe the phenomenology of the coexis- 
tence of d-wave superconducting and antiferromagnetic 
regions. This approach is similar to the starting point of 
many recent Bogoliubov-de Gennes calculations^^^^^^^^. 
The Hamiltonian in Eqn. (l)-(2) can be viewed as 
the mean-field Hamiltonian oi a t — U — V Hubbard 
model, where the nearest neighbor attraction V gives 
rise to the d-wave superconductivity. In contrast the on- 
site Coulomb repulsion U only causes the antiferromag- 
netism. In this article we do not diagonalize H in the 
Bogoliubov-de Gennes scheme since such lattice calcu- 
lations require unrcalistically large gap A and magnetic 
field values. Instead we solve the Dyson equation exactly 
by inverting a large matrix. This approach has previ- 
ously been utilized extensively to study various short- 
ranged impurity effects in superconductors'^'^, but can 
also be used for extended perturbations embedded in a 
Go medium. Here Go is the Green's function of the par- 
ent medium, in this case a d-wave BCS superconductor. 
This Green's function is given by 

Go ^ (p, w) = (w + i5)To - ^pT3 - ApTi (3) 

where denote the Pauli matrices in Nambu space and 
the gap function Ap = ^ (cos {px) — cos {py)). The lat- 
tice constant oq is set to unity and ^p = Cp — M with 

gp = - 2t (cos {px ) + cos (py))-U' (cos ) cos (py ) ) . (4) 

Here t{t') refers to the nearest (next-nearest) neighbor 
hopping integral and is the chemical potential. We 
perform the 2D Fourier transform of Go(p,cj) numeri- 
cally by utilizing a real space lattice of 1000 x 1000 sites 
and a quasiparticle energy broadning oi 5 = l.OmeV. 
To simulate the situation around optimal doping of the 
hole doped cuprates the following parameters are cho- 
sen: t = 300meV, t' = -120meV, Aq = 25meV, 
/U = — 354meV. When the real space domain affected 
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by iJ™* involves a finite number of lattice sites N x N 
we can solve the Dyson equation exactly to find the full 
Greens function. Writing the Dyson equation in terms of 
real-space (and Nambu) matrices it becomes 



(5) 



The size of the matrix (l - ii*"*G°(ti')) is given by 
{d X N^) X {d X N"^) where d is an integer equal to the 
number of components in the Nambu particle-hole spinor 
and N denotes the total number of lattice sites affected 
by the magnetic perturbation. Therefore a real-space lat- 
tice with 25 X 25 sites affected by perturbations results 
in a (1250 x 1250) matrix to be inverted. 
Knowing the full Greens function we obtain the LDOS 
p{r,LLi) = -ilm [Gii(r, w) -I- G22(r, -w)], which is pro- 
portional to the differential conductance measured in the 
STM experiments. 

Wc have checked that the above approach reproduces the 
expected LDOS for unitary non-magnetic impurities in d- 
wave superconductors^". Also in this one- impurity case 
we reproduce the constant-energy LDOS maps recently 
calculated by Wang ct al.^^'^"^ 

Motivated by the qualitative agreement of the spin struc- 
ture in figure 2a with the neutron data, we assume that 
this represents the induced magnetism around the vor- 
tices and calculate the LDOS in this striped environment. 
To this end we simply restrict the sum in Eqn. (2) to 
include the sites within these magnetic regions. The sys- 
tem is depicted in figure 2a where the grey background 
reveals the underlying superconducting state. Again the 
black (white) squares correspond to the sites affected by 
the staggered magnetic perturbation. 
Figures 3 and 4 show real-space maps of the LDOS 
summed over a small energy window from -12 meV to 
+12 meV in intervals of ImeV for different strengths of 
the antiferromagnetic perturbation M. The vortex cen- 
ter is located in the center of the images. Figure 3 (4) 
is calculated with (without) the antiphase modulation of 
the adjacent stripes. Thus the spin configuration of fig- 
ure 2a corresponds to the images in figure 3. The clear 
difference between the LDOS images of figures 3 and 4 
reveals that the STM technique can be used to determine 
this phase relation. It is clearly seen from both figures 
3 and 4 that the low energy LDOS structure eventually 
becomes ringshaped as the magnitude of M increases. 
In this limit the pinned stripes operate as steep poten- 
tial walls. Figures 3a and 3b seem to display the closest 
reseniblence to the experimental data^^ which indicates 
that the induced magnetism is very weak. In figure 5 we 
show the Fourier transform of several constant energy 
LDOS images for M = 100 meV with the antiphase spin 
modulation included. In these figures the Fourier compo- 
nent q = is located at the center. The detailed energy 
dependence of these images is caused by quasiparticlc in- 
terference effects as pointed out by Wang at al.^^ in the 
case of a single impurity. 
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FIG. 3. Real-space LDOS summed from -12 meV to +12 

meV for: a) M = 35 meV, b) M = 100 meV, c) M = 200 meV, 
d) M = 300 meV. 



The dispersive features of the images presented in figure 
5 are dependent on the microscopic parameters and the 
associated Fermi surface. However, it is also evident that 
the ringshaped stripes surrounding the vortex cores give 
rise to non-dispersive intensity around q = —(±1/4,0) 
and q = ^(0, ±1/4). This in turn leads to the checker- 
board pattern in the low energy sums of the LDOS dis- 
played in figures 3 and 4 whereas the dispersive features 
fade away in these summed LDOS images^^. 




FIG. 4. Rcal-spacc LDOS summed from -12 moV to +12 
meV for: a) M = 35 meV, b) M = 100 meV, c) M = 200 meV, 
d) M = 300 meV. 
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We have confirmed this fact by identifying similar non- 
dispersive features in the LDOS around configurations 
with different periodicities. For instance a structure with 
2ao charge periodicity leads to non-dispersive intensity 
around q = f^(±l/2, 0) and q = ^(0, ±1/2). 




non-dispersive feature in the LDOS arising from the in- 
duced static antiferromagnetism. This feature gives rise 
to the checkerboard LDOS observed experimentally by 
Hoffman et al}^ Thus both the STM measurements and 
the enhanced intensity of the incommensurable peaks ob- 
served by neutron diffraction can be ascribed to the phe- 
nomena of a single Cu02 plane. 

Acknowledgement: Support by the Danish Natural Sci- 
ence Research Council, Ole R0mer Grant No. 9600548. 




FIG. 5. Fourier images of the constant energy LDOS maps 
for M — 100 meV and a) a; = 3 meV, b) a; = 6 meV, c) a; = 
9 meV, d) w = 12 meV. 

In the above calculation wc have not yet included the 
Doppler shift from the circulating supercurrents or the 
gap suppression close to the vortex core. As pointed out 
by Polkovnikov ct al}^ the former effect is not expected 
to produce significant changes of the four-period mod- 
ulations. As for the latter we have checked that a gap 
suppression only leads to minor quantitative changes in 
the dispersive part of the LDOS. Finally, Podolsky et 
ai."^^ discussed scenarios of weak translational symmetry 
breaking and found that in order to explain quantita- 
tively the zero-field STM results by Howald et aJ.^^ one 
needs to include dimerization, the modulation of the elec- 
tron hopping. This dimerization will also produce quan- 
titative changes, but not alter the qualitative conclusion 
that pinned stripes produce checkerboard LDOS. 
In summary we have discussed the phenomenology of a 
simple physical picture of pinned stripes around vortex 
cores which are forced to be antiferromagnetic by an ap- 
plied magnetic field. The induction of magnetic striped 
race-tracks around the core is consistent with the neu- 
tron diffraction spectra observed on LSCO with a dop- 
ing level near x=0.10. As expected this is only true if 
the stripes are out of phase with their neighbors in the 
usual sense. In materials where a checkerboard spin pat- 
tern is relevant (possibly Bi2212 or overdoped LSCO), we 
show that a 45 degree rotation of the main incommen- 
surable peaks is to be expected. Finally we studied the 
electronic structure around the vortices and identified a 
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